THE SPECTRUM OF THE HILBERT SPACE VALUED SECOND DERIVATIVE 
WITH GENERAL SELF-ADJOINT BOUNDARY CONDITIONS 



JOACHIM VON BELOW AND DELIO MUGNOLO 

Abstract. We consider a large class of self-adjoint elliptic problem associated with the second derivative 
acting on a space of vector-valued functions. We present two different approaches to the study of the associated 
eigenvalues problems. The first, more general one allows to replace a secular equation (which is well-known 
in some special cases) by an abstract rank condition. The latter seems to apply particularly well to a 
specific boundary condition, sometimes dubbed "anti-Kirchhoff" in the literature, that arise in the theory of 
differential operators on graphs; it also permits to discuss interesting and more direct connections between 
the spectrum of the differential operator and some graph theoretical quantities. In cither case our results 
yield, among other, some results on the symmetry of the spectrum. 



1. Introduction 

Differential equations on networks have a long history, starting probably in 1847 with Kirchhoff 's electrical 
circuit equations using the potential mesh rule and an incident currency law [3T]. For specific differential 
edge operators both latter conditions were naturally replaced were by the continuity potential condition and 
an incident flow condition. Nowadays, dealing with rather general vertex transition condition, we are led to 
abstract interacting problems, as e.g. the Bochner space setting below. 

The problem of determining all possible self-adjoint realizations of a differential operator is quite common 
in mathematical physics, as it is always possible conceivably to prepare a corresponding quantum system in 
such a way that the associated Hamiltonian is one such realization. 

Relying upon known results for differential operators on domains, in [22 V. Kostrykin and R. Schradcr 
have proposed a natural representation of self-adjoint boundary conditions for 1-dimensional systems, i.e., 
for Laplace operators on a Bochner space L 2 (0, 1;£ (E)). If E is the edge set of a graph, then such systems 
are usually called "networks" or "quantum graphs" in the literature. That approach has the drawbacks 
that it is poorly fitted for a variational setting, and that the boundary conditions do not determine their 
representations uniquely. Both issues can be avoided making use of an alternative parametrization proposed 
by P. Kuchment in [23] ■ Both parametrization are equivalent, in the sense that each Kuchment's boundary 
condition can be represented using Kostrykin-Schrader's formalism, and vice versa. Moreover, if E is finite, 
i.e., if one considers a differential operator on finitely many edges, then it has been shown in [53] that self- 
adjoint boundary conditions are actually exhausted by either parametrization. Kuchment's parametrization 
can be reduced to the choice of a closed subspace Y of £ 2 (E) and of a bounded linear self-adjoint R on £ 2 (E); 
see Section[2]for details. To the best of our knowledge, not much is known about spectral properties of general 
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sclf-adjoint realizations, beside the resolvent formula of Krein type obtained in [28 | 129 ] . The goal of our note 
is to determine the (real) spectrum of this class of Hamiltonians in dependence of Y and R. 

It turns out that, on a spectral level, a different choice of R has the sole effect of shifting (in a complicated 
way) the eigenvalues, but the essential choice is, in a certain sense, that of Y. This becomes particularly 
clear at the end of Section [2] where, after obtaining a general description of the spectrum in the general case, 
we present a particularly appealing spectral symmetry if R = 0. Moreover, we show how our results can be 
further simplified in the case of small sets E - and in particular in the case of £ 2 (E) = C,i.e., of a single loop. 

There are at least two special cases where our abstract secular equation has a much more intuitive meaning: 
the boundary conditions may reflect so-called continuity and Kirchhoff conditions - this is in fact by far the 
most common condition considered in the literature, cf. Section [3] In Section [4] we consider a further - less 
standard but interesting - boundary condition, which is in some sense dual to the previous one: the function 
satisfies Kirchhoff 's law where its associated flow is continuous in the network's vertices. 

As a motivation for the study of this special and seemingly counter-intuitive boundary condition, it should 
be remarked that this so-called "anti-Kirchhoff" condition is already known to be distinguished among all 
possible self-adjoint conditions for the second derivative on a graph. E.g., it has been observed in [T^l [37] that, 
apart from the decoupled conditions of Dirichlet/Neumann type, Kirchhoff and anti-Kirchhoff conditions are 
the sole ones invariant under edge permutations (up to lower order perturbations); and on a loop, also the 
sole ones that induce L°°-contractive heat semigroups. 

In Sections [3] and [4] we study these two boundary conditions thoroughly. In this special case, the general 
secular equation turns out to be tightly related to the characteristic polynomial of the transition matrix of the 
underlying abstract graph. We can then strongly improve our general results and are in fact able to describe 
the spectrum of the Laplacian under said conditions in terms of invariants of the underlying graph: several 
features are proved to reflect directly known graph theoretical properties and quantities, like bipartiteness. 
Our results can be seen as partial solutions to certain inverse problems. In Section [5] we compare information 
made available by our results with the outcome of previous linear algebraic investigations. 



2. Eigenvalue problems for the second derivative on general Bochner spaces 

Let H be a separable Hilbert space and consider the Bochner space L 2 (0, 1; H). Let Y be a closed subspace 
of H x H and consider the Hclmholtz equation with boundary conditions 



(2.1) 

i.e., the eigenvalue problem 



(EP Y ) 




x e (0,1) 



in L 2 (0, 1; H). This kind of boundary conditions might be generalized considering 



(EPy.r) 



R 




G 
G 



Xu(x), 



x G (0,1) 



Y A 



«'(i) j 1 " wi), 

where R £ £{¥)■ This class of boundary conditions has been made popular by [23], but we point out that 
they, and even some generalization to eigenvalue-dependent boundary conditions, had already been suggested 
in |1 §3]. 
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Remark 2.1. In the special setting of diffusion on I -dimensional ramified structures (so-called networks 
or quantum graphs,) this class of boundary conditions has been proposed in |23j . Observe that the boundary 
conditions in (EPy,r) are less general than those of Agmon-Douglis-Nirenberg type 

(ADN) AftlV^flS^O 



u(l)J " V 

for A, B G C(H x H): this is due to the fact that the first boundary condition in (EPy.r) determines the 
second completely. E.g., if we take H = C, then taking 

i.e., letting 

/(0) = /(l) 

forces the second boundary condition to read 

/'(0) = /'(l) + i?/(0) 

for some R £ C . Thus, if we consider 

A:=(l "M and B := ( ° _° 



thus obtaining continuity 

and a second boundary condition 



/(0) = /(I) 



2/'(0) - /'(I) 



other than Kirchhoff, these both boundary conditions cannot be expressed in the form of (EPy,r) for any 
choice of R. 

In spite of their less general character, the boundary conditions in (EPy,r) give the problem an interest- 
ing variational structure that will prove extremely useful in spectral investigations. The sesquilinear form 
associated with (EPy,r) is given by 

<**(/, 3):= [\f(x)\g\x)) H dx+(R^ m ^ ^ (0) 



with dense domain 



Vr:={feH\o,i;H)\ e y} . 

We denote by —Ay,R the operator associated with (or, Vy), i.e., 

D(A Y . R ) := {feV Y \3geL 2 (0,l;H)s.t.a R {f,h) = (g\h) H iov a \\heV Y } 

- {fe B^,1;B)\ (-/$>) + « Q e^«d (per 
Ay.ii/ := -3- 

The model is still sufhciently easy to allow us to explicitly compute the eigenvalues of the corresponding 
problem (EPy,r). We mention that, with a different representation of the boundary conditions, the number 
of negative eigenvalues has been computed in [7J. 
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Proposition 2.2. Let Y be a closed subspace of H x H and R € C(Y) be self-adjoint. The spectrum of the 
operator —Ayr is real; it is contained in [0, oo) if R is positive semide finite. If H is finite- dimensional, then 
— Ay,R has pure point spectrum. Then the following assertions hold. 

1) is an eigenvalue of (EPy,r) if and only if 

H Y :={(A,B)eHxH:A = B}r\Y^{0} and R = 0, 

and in this case the multiplicity of agrees with dim Hy < dim H. 

2) Let A > 0. If R is positive semidefinite, then A > is an eigenvalue of (EPy,r) if and only if the space 
Hym of all solutions (A, B) e H x H of the system 



(2.2) 



Py± (acos\/A + Bs'mVx) ° 

Py f VX f J 4 sinx /T_ 5 COSA /X J R v Acosv / A + B sin V\ 



has nonzero dimension; and in this case the multiplicity of X agrees with dimiJy jjR < 2dimiJ. 

Proof. The fact that the form (an, Vy) is real and, for positive semidefinite R, accretive shows that the 
spectrum of — Ay r is contained in the positive half line. If H is finite-dimensional, then the compact 
embedding of D{Ay^r) into L 2 (0, 1; H) follows from the Aubin-Lions Lemma and yields that — Ay^R has pure 
point spectrum. 
(l)If 

u"(x) = o, xe[o,i], 

then 

o = a {u ,u) = jf 1 \\A*)w%d*+ (r (:[?|) i ($)) HXB - 1! 

hence u can only be an eigenvector for the eigenvalue if it is entry-wise constant, i.e., 

u(x) = AeH for all x e [0,1]. 

Hence, in order for u to satisfy the first boundary condition in (EPy,r) (the second one being satisfied if and 
only if R = or Y = {0}, of course) one needs that 



«(0)\ = A 
u(l) - [A 



G Y. 



(2) For A > 0, we see that the characteristic equation 

-u"{x) = Xu(x), x e [0, 1], 

is solved for all A > by 

u(x) := A cos VXx + Bs'm VXx, i£ [0,1], 

for some A, B £ H, whence 

u'(x) = VX(— Asin^/Xx + Bcos\^Xx), x e [0, 1], 

and hence 

'm(0)\, _( A 



along with 

-u>(0)\ -B 



u(l)J \A cos y/X + B sin VX 

0) ) =Sx( - B 

u'(l) I \-Asin^X + BcosVX 
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Imposing the boundary conditions (2.1 1 we obtain 

A 



and 



A cos v A + B sin v A 



Py (^(ylsinVA-BcosV^) ^ cos %/A + B sin y/\J ) °' 

Hence, A > is an eigenvalue if and only if the above equations have a nontrivial solution (A, B) £ H x H. 

Combining the cases in (1) and (2) yields the claimed rank condition. Finally, the bound on the dimension 
of the eigenspace associated with follows from Grassmann's formula, as 

dim H Y = dim H + dim Y - dim({(A, B) e H x H | A = B} + Y) < dim H. 

This concludes the proof. □ 

If R = 0, an easy but interesting consequence of this general result is the following. 

Corollary 2.3. Let Y be a closed subspace of H x H and R — 0. Let both A > and (tt — y/X) 2 > 0. Then A 
is an eigenvalue of (EPy) if and only if (tt — V~X) 2 is an eigenvalue of (EP Y i). In this case, their algebraic 
multiplicities coincide. 



Proof. By Proposition 2.2 A is an eigenvalue of (EPy) if and only if the space Hy of all solutions (A, B) 6 
H x H of the system 

Y± I Acos v / A + Bsinv 7 ^/ ~° 



(2.3) 

has dimension larger than 0. Now, it suffices to observe that because 



A sin y/X — B cos Vx - ' 



sin y (tt — V~X) 2 = sin vA and cos y (n — VX) 2 = — cos y/X 



(A,B) is a solution of the system in (2.3) if and only if (A,B) := (B,A) is a solution of the system in (2.3) 
with the roles of Y, Y 1 - interchanged - i.e., of 



(2.4) 



P ( 1=0 

\A cos vA + B sin yA j 

Y \A sin y/X — B cos y/X , 



This concludes the proof. □ 

Example 2.4. In the special case where (EPy.r) can be represented as an eigenvalue problem on a network 
associated with a directed graph T, the most natural boundary conditions are those requiring continuity of the 
function values as well as a Kir chhoff -like rule for the normal derivatives at each vertex ofT, see Section^for 
details. Such conditions have been discussed in the literature since [25] , and in fact implicitly even since |31) . 
It is well-known that the associated elliptic problem is self-adjoint, and in fact we can represent them in our 
formalism letting 

Y := Range I?, 

where 

(v-f 



(2-5) •= I it)—\t 



6 



JOACHIM VON BELOW AND DELIO MUGNOLO 



and T> + ,T>~ are the matrices whose entries are the positive and negative parts of the entries of the (signed) 
incidence matrixT) ofT . In this case one can check that the dimension ofY agrees with the number of vertices, 
i.e., with the cardinality of V , but for general Y there is no such an intuitive interpretation. An exception is 
constituted by the "dual" conditions to the above one, i.e., by the conditions arising from the choice 

Y := (Range V) ± , 



which in the literature are usually referred to as (5'-couplings or anti-Kirchhoff vertex conditions. 



Definition 2.5. For V defined as in (2.5) and with 



Y := Range V 



Y := 



(Range Vj 



the boundary conditions 



Y. 



A 1 ) < 
-u'(oy 

u'(l) 

will be referred to as CK-condition and KC -condition, respectively. The second derivative with these condi- 
tions will be denoted by A CK and A KC , respectively. 

Several features of the JiTC-condition - including index theorems, trace formulae, convergence results, 
parabolic properties - have been considered among others in [HI [TJ [2Ql |30l [13] . We will thoroughly discuss 
the eigenvalue multiplicities for these boundary conditions in Sections [3] and [4] 

Remark 2.6. If H is finite-dimensional, say dim H = m, and R is positive semidefinite, then the Am x 2m 

r , ( ) \ 

Y \cos y/X Idem sin y/X Idc»» / 

/ n TJ . \ / J A _. Pi \ \ 

Py Vx 



algebraic system (2.2) has a nontrivial solution (A, B) if and only if 



( 



(2.6) 2m — rank 



V 



o 



sm 



Id cm 
COS y/X Idc 

Id cm 

cos vA Idc™ 



sin v A Idc™ / 

Idem 

^COS 





sin y/\ '. 



> 1. 



lim , 



Hence, by Proposition 2.2 each eigenvalue X can have multiplicity m(X) at most 2m. More precise bounds on 
m(A) can be obtained in the special cases considered in Sections [5] anrfj^J 

As the eigenvalues A^ interlace with the Neumann and Dirichlet eigenvalues on N uncoupled intervals a). 
and Wfc respectively, i.e. ak < Afc < ^>k for all k € N, we can state the 

Corollary 2.7. For the eigenvalues of the Laplacian associated to the problem (EPy,o) the following spectral 
asymptotic hold 

A* _ ^ 
N 2 ' 

where X^ denotes the k-th eigenvalue, and where the eigenvalues are counted according to their (geometric = 
algebraic) multiplicities. 

Let now R = 0. We close this section focusing on the special case of H = C, i.e., of an elliptic problem 
concerning only one interval. As long as we want to keep locality of the boundary conditions, there are only 
four possibilities: 

• Dirichlet boundary conditions in both and 1; 

• Neumann boundary conditions in both and 1; 

• Dirichlet boundary condition in and Neumann boundary condition in 1; 

• Neumann boundary condition in and Dirichlet boundary condition in 1. 
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They correspond to the boundary conditions in (EPy) by means of four different subspaces Y of H x H = C 2 : 
these are respectively 

• Y = {0}; 

• Y = C 2 ; 

• Y = subspace spanned by the vector (0, 1); 

• Y = subspace spanned by the vector (1,0). 

If we generalize the above setting in order to allow for nonlocal interactions between the vertices at and 1, 
we are actually performing an identification between two vertices: we might as well regard this setting as a 
differential operator on a loop around a vertex v. 

The boundary conditions in (EPy) are now given by a subspace of C 2 : besides the same four boundary 
conditions appearing in the interval case there are infinitely many new ones. All such boundary conditions in- 
terpolate between the Dirichlet-Dirichlet and the Neumann-Neumann case, i.e., Vy contains the form domain 
corresponding to the the Dirichlct case (i.e., ffg(0, 1)) and is contained in the form domain corresponding 
to the Neumann case (i.e., 7? 1 (0, 1)). In particular, the operator associated with (ao,Vy) is dominated by 
(resp., dominates) the second derivative with Dirichlet (resp., Neumann) boundary conditions, in the sense 
of self-adjoint operators. 

It seems that only a few publications are devoted to the study of differential operators on a loop. Among 
them, we mention |24j . where a characterization of a certain class of second order self-adjoint elliptic operators 
is presented. 



Already Proposition 2.2 shows that the second derivative Ay,o with associated eigenvalue problem (EPy) 
is not invertible if and only if Y contains vectors of the form (A, A) for some A 6 H. In the case H = C 
this means that Ayo is not invertible in exactly two cases: Y = C 2 and the space Y spanned by the vector 
(1,1), corresponding to (uncoupled) Neumann boundary conditions continuity/Kirchhoff conditions at the 
endpoints, respectively. In this section we are going to describe the spectrum of (EPy) more explicitly. 

We focus on the boundary conditions defined by a 1-dimensional subspace Y and thus neglect the trivial 
cases of Y = {0} 2 or Y = C 2 (standard Dirichlet or Neumann boundary conditions, whose associated spectrum 
is well-known). We can then consider conditions of the same form as in (EPy) with 

Y = Y a := / (f\ \ and Y ± = Y^ = Y_ & -x 



The associated orthogonal projections are given by 



(2.7) P Y = k ( 1 i J and P Y ± = k . ( 



a IN , „ /4 -1 



for 



Re a 

k := 



1 + Re a 

Remark 2.8. In order to consider the totality of possible boundary conditions it seems at first that we have 
to allow all a G C U {oo}, but in fact the eigenvalue problem (EPy Q ), for some \a\ > 1, can be equivalently 
realized (i.e., we have isospectrality) by applying the parity transformation 

u(l/2 + x) h-> u(l/2 - z), a; 6 [0,1], 

and then switching to the eigenvalue problem (EPy _ 1 ). 

Proposition 2.9. For all a € C the spectrum of (EPy Q ) is given by 

, „. ^ 2Rea 
A>0|cos^ =rTRF 

All the eigenvalues are simple for all a ^ ±1. All the eigenvalues have multiplicity 2 for a — ±1. 
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Proof. First of all let us observe that, by Proposition |2.2[ is an eigenvalue if and only if a = 1 



Now, plugging (2.7) into (2.6) and eliminating the linearly dependent rows yields that A > is an eigenvalue 



if and only if the matrix 
(2.8) 

is singular, i.e., if and only if 



— a cos —a sin vT 
sin yX a — cos y/X 

2 Re a = cos\/A(l + |a| 2 ). 



Moreover, the eigenvalue A has multiplicity 2 if and only if the matrix in (2.8 1 vanishes: this is the case if and 



only if cos vA = 1 and a = 1 , or else cos vA = — 1 and a = — 1 . □ 



Remark 2.10. It follows from Proposition 2.9 that the spectrum associated with Y\ (CK-condition, with the 
notation introduced in the Example 2.4) is given by {(2fc7r) 2 | k € No} for a — 1, while the spectrum associated 
with Y_x (KC- condition) agrees with {((2k + l)ir) 2 \k G No} for a = —1. This is of course nothing but a 
special case of Corollary \2.S\ Hence, among all second derivatives with boundary conditions as in (EPy q ), 
those corresponding to a = 1 (i.e., to the CK-condition) have the smallest possible lowest eigenvalue; and 
therefore, those corresponding to a = —1 (i.e., to the KC -condition) have the largest possible lowest eigenvalue. 
Of course, these considerations have a counterpart in the asymptotics of the associated semigroups that govern 
the heat equation corresponding with (EPy Q ). 

3. Adjacency calculus and spectral analysis under the CRT-condition 



In this section we recall some known results [2j [4] on the Laplacian on a graph by means of combinatorial 
quantities when continuity and Kirchhoff type conditions are imposed in the vertices. For that purpose, 
let us introduce some terminology. For any graph T — (V,E,E), the vertex set is denoted by V = V(T), 
the edge set by E = E(T) and the incidence relation by ec V x E. The degree of a vertex v is defined 
by j(v) = card{e £ E\v € e}. All graphs considered in this and the following section are assumed to be 
nonempty and finite with 

n = #V, N = #E 
and c = c(r) connected components Ti, . . . , T c . 

Definition 3.1. 

c + (r) = #{r fe |r fc bipartite} 
c~(T) = f/:{Tk\Tk non-bipartite} 

The graphs T are also assumed to be simple, i.e. T contains no loops, and at most one edge can join 
two vertices in T. By definition, a circuit is a connected and regular graph of degree 2. We number the 
vertices by Vi, . . . , v n , the respective degrees by 71, . . . , 7„, and the edges by ei, . . . , ejy- The adjacency matrix 
A(T) = (e,/j) of the graph is defined by 



1 if Vi and Vh are adjacent in V 
else 



Note that *4(r) is indecomposable if and only if T is connected. By simplicity, any two adjacent vertices Vi 
and Vh determine uniquely the edge e s joining them, and we can set 



s if e s C\V = {vi,v h }, 
1 otherwise. 



s(i, h) 

For further graph theoretical terminology we refer to [31] , and for the algebraic graph theory to [HI US] ■ 



SPECTRAL ANALYSIS OF THE SECOND DERIVATIVE 



9 



The abstract graph T can be realized as a topological graph in K m with m > 3, i.e. V(T) C K m and the edge 
set consists in a collection of Jordan curves E(T) — {ttj : [0, 1] — ^ M| 1 < j < N}. The trace of the topological 
graph leads to the associated network or metric graph G — UjLi n j (IP> !])■ The arc length parameter of the 
edge ej is denoted by Xj. Clearly, for many purposes the irj have to display certain regularity properties, 
say Wj £ C 2 , but for the present context it suffices to identify each edge ej with [0, 1]. We shall distinguish 
the boundary vertices Vb = {vi € V\ji = 1} from the ramification vertices V r = {vi € V|7» > 2}. The 
orientation of the graph T is encoded in the signed incidence matrix, which in this setting can be written as 

if 7Tj(l) = Vi, 

dij = ^ — 1 if 7Tj(0) = Vi, 
otherwise. 



In fact, the factors dij stand for the outer normal derivative at Vi with respect to ej. Functions on the graph 



or on the network are mappings u = (uj) Nxl ■ [0, 1] — > C N with edge components Uj : [0, 1] — > C. Here we 



and use the abbreviations 

. .. a . . 

etc. 



Q 

Uj{v t ) := Uj(nJ (vi)), djUj(vi) := —Uj{xj) 



It is well-known that corank(r) := dime ker 2?(r) = N — n + c(T), and that corank(r) = 1 if and only if T is 
unicyclic, i.e. by definition, T contains exactly one circuit. In the connected case, T is unicyclic if and only if 
N = n. 

In this short section we consider the standard Laplacian 



v 



A = ((uj) Nxl ^ (%Uj) Nxl ) : C 2 [0, 1} N C[0, 1] 
under the continuity condition 

(3.1) Wi € V r : ej n e s = {v,} Uj(vi) = u s (vi), 
and under the Kirchhoff flow condition 

N 

(3.2) '^^d i jdjUj(v i ) — for 1 < i < n. 

As anticipated in Section [2j we shall refer to Conditions (3.1) and (3.2 1 as to the C K -condition. It can be 
readily written as a canonical boundary condition in the sense of Holder as follows, see [3]- Define the matrix 



-Poo Pqi 
Pio Pu 



2JVx2Af 



with matrices P a p — (p a /3jk)NxN defined by 

Pafijk = 



1 if TTj(alj) = n k (pi k ), 
otherwise. 



The matrix S2 is symmetric and decomposable into a block diagonal matrix with n blocks, each of them being 
a dyad of the form e 7i e* , where throughout = (l)fcxi- Let Si denote the orthogonal projection onto the 
kernel of S2. In detail, for dij = —1 or dij—N = +1 with 1 < j < 2N and z — (zj)2Nxi we have 

(Siz)j = Zj - — ^2{z k \ 1 < k < 2N,d lk = -1 or d itk _ N = +1} . 
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Xow I bo continuity condition at ramification vort ices roads Si ( | = I), and the ( ' K Condition lakes tin- 
form 



ker S, 



^(i)J k 1; U'(l) 

This shows that we are in fact inside the general setting presented in Section [2j It is well known that the 
eigenvalues of —A are nonnegative, see e.g. Proposition 2.2 or [2 . Following the transformations in [213] the 
eigenvalue problem for A is equivalent to the matrix differential boundary eigenvalue problem ( |3.5[ )-( |3~T0 ) 
below incorporating the adjacency structure of the network. For that purpose we recall that the Hadamard 
product of matrices of the same size is defined as (a^) * (bik) = {aikbik)- For a function u : G — >■ R denote its 
value distribution in the vertices by 



(3.3) 

and for x G [0, 1] define the matrix 



ip = n(u) = 00i)) nxl 
V{x) = {u ih {x)) nxn 



by 

(3.4) 



u ih (x) = e ih u s (i^ h ) 



s(i,h) 



— xd 



is(i,h) 



Then the eigenvalue problem for A under (3.1) and (3.2) reads as follows 



(3.5) 


u ih eC 2 ([0,1]) 


for all 1 < i, h < n 


(3.6) 


eih = =>■ u ih = 


for all 1 < i, h < n 


(3.7) 


U" = -AU 


in [0, 1] 


(3.8) 


U(0) = ip e* *A 


continuity in V^.(r) 


(3.9) 


U*(.x) =U(l-s) 


for x € [0, 1] 


(3.10) 


U'(0)e = 


Kirchhoff flow condition 


Set 







$:=TJ(0) = <pe**A, *:=U'(0). 
Finally, introduce the row-stochastic transition matrix 

(3.11) Z := Diag(^ e)" 1 (A), 

that has only real eigenvalues [2] and plays a key role in the characterization of the spectrum of the Laplacian. 
As for the multiplicities recall the following, which has been obtained in [2JH]. 

Theorem 3.2. If X is an eigenvalue of —A inC\{G) and ip € K™ is a vertex distribution of an eigenfunction 
belonging to \, then 

(3.12) Zip — cos vA tp. 

Conversely, if cos y/X is an eigenvalue of Z admitting the eigenvector ip € K" 7 then X is an eigenvalue of —A 
inC\(G) and ip the vertex distribution of some eigenfunction belonging to X. The multiplicities are 



m g (X) = m a (X) 



'c(r) 

m ff (cos vCa, Z) — m a (cos vA, Z) 
corank(r) + c(T) = N -n + 2c{T) 
,corank(r) + c+(T) - c"(r) = N - n + 2c+{T) 



ifX = 0, 

if sin vA 7^ 0, 

if cos vA =1, A > 0, 

if cos VA = — 1, A > 0. 
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Remark 3.3. In other words, in the relevant case of a connected graph the spectrum of —A under the CK- 
condition can be partitioned as 

<t(-A) = {0} U a, U a si 

where 

o~i := |(2^7r ± arccosa) 2 | a € cr(Z) \ {— 1, 1} and I 6 Z j , 

and 

a s ■= {k 2 n 2 \k e Z\ {0}}. 

As the elements of Ci are solely determined by the adjacency structure of the graph, they are sometimes called 
"immanent eigenvalues". For the multiplicities of the eigenvalues we have for all k G Z \ {0} 

(1) m(0) = 1; 

(2) m(A) = dimker ^cos \f\ + Z^j for A £ c^; 

(3) m(k 2 ir 2 ) = N — n + 2 ifT is bipartite and m((2fc + l) 2 7r 2 ) = N — n ifT is not bipartite; 

(4) m(4fc 2 vr 2 ) = N - n + 2. 

A spectral asymptotics of Weyl type follows promptly. 
Corollary 3.4. ('[HH],) For the canonical Laplacian —A under the CK -condition the spectral asymptotics 

y \ k 7T 2 

holds, where Afe denotes the k-th eigenvalue, and where the eigenvalues are counted according to their (geo- 
metric — algebraic) multiplicities. 



By Theorem 5.6 below, a necessary conditions for the whole spectra of — A CK for two different graphs 
Ti and T2 to be comparable is that both graphs have the same number of vertices, edges, bipartite and 
non-bipartite components. This is of course a very strong assumption. Thus, a more natural question is 
whether some relevant spectral subset - i.e., the lowest non-trivial eigenvalue - of —A CK is lowered or raised 
by certain graph operations. For instance, the following comparison result for eigenvalues of A CK is a direct 



consequence of Theorems 3.2 and 4.1 (4.161 and a result by Chung jTSJ Lemma 1.15 



Proposition 3.5. Let Ti,r2 be two connected graphs such that T2 be formed by contractions of vertices from 
T\. Then, the lowest nontrivial eigenvalue of —A CK on the metric graph G\ associated with Y\ is at most as 
large as the lowest nontrivial eigenvalue on the metric graph Gi associated with T2 ■ Hence, the heat equation 
under CK -conditions on G\ converges to the equilibrium faster than on G 2 - 

4. Spectral analysis of the ifC-coNDiTioN 
In this section, we consider the standard Laplacian 

A = (K-) Wx i ^ (^).vxi) : C2 t°' l f -+ C ^ ^ 
under the so-called anti-Kirchhoff condition, i.e. 

N 

(4.1) d ij u j( v *) = for 1 < * < n 

3=1 

and 

(4.2) Vwi € V r : Bj He s = {uj =^ d^u^Vi) = d is u s (vi), 
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As in Definition 2.5 we refer to (4.1) and (4.2) as the KC -condition. In terms of the conditions discussed in 



Section [2j they correspond to 



u(0) 
u(l) 



G Y 



kerS, 



-«'(0) 
u'(l) 



ker Si, 



where Si,S2 are defined as in Section [3j It follows from Proposition |2.2| that the eigenvalues of —A under 
these transition conditions are nonnegative. Using the adjacency transforms from Section |3j the eigenvalue 
problem for A reads as follows. 



(4.3) 
(4.4) 
(4.5) 
(4.6) 
(4.7) 
(4.8) 



u%h eC 2 ([0,i}) 

e lh = 0=> u ih = 
U" = -AU 
U*(a;) =U(l-a;) 

U(0)e = 
U'(0) = ^e**^ 



for all 1 < i, h < n 
for all 1 < i, h < n 
in [0, 1] 
for x G [0, 1] 



(4.1) 



(4.2) 



As above, set $ := U(0), * := U'(0) = ip e* *A and Z = Diag (A e) A. Then we can state the following. 

Theorem 4.1. 1) If A > is an eigenvalue of —A under the KC-condition in C 2 [0, 1} N and ijj G M n is a 
vertex distribution of the normal derivatives of an eigenfunction belonging to \, then 



(4.9) 



Conversely, if A > and — cos yA is an eigenvalue of Z admitting the eigenvector ip G M™, then A is an 
eigenvalue of —A under the KC-condition in C 2 [0, 1] N and ip the distribution of normal derivatives of some 
eigenfunction belonging to A. 

2) The multiplicities of the eigenvalues are given by 



m g (X) 



'corank(r) - c"(r) = N - n + c+(T) 
i g (— cos vA, Z) = m a (— cos vA, Z) 
corank(r) + c(T) = N-n + 2c(T) 
^corank(r) + c+(T) — c~(T) = N — n- 



2c+(T) 



ifX = 0, 
if sin 

if cos VA = —1, A > 0, 
if cosy/X = 1, A > 0. 



Let us recall from [2] that for T connected, i.e. c(T) = 1, 

dim R 7W _ (r) = corank(r) = N — n + 1, 

' N - n + 1 if T bipartite, 
N — n if r is not bipartite, 



(4.10) 
(4.11) 
where 



dim R X + (r) 



M(T) {M\M = {m ih ) nxn , Vi, h G {1, . . . , n} : e lh 
M~{Y) := {M G M{Y)\M* = —M, Me = 0} , 
M + (T) := {M £ M(T) \ M* = M, Me = 0} . 



m lh = 0} , 



Proof of Theorem \4-.l\ As corank(r) = J2k=i corank(Ffe), we can confine ourselves to the case c(T) = 1. First, 
let us consider the case 

A = 0. 
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Then each eigenfunction of (4.3 1-( |4.8[ ) is of the form 

(4.12) U(x) = $ + a;* with * £ M(T). 

Reasoning as in the proof of Proposition |2.2[ the slope matrix has to vanish, since u has to be constant on 
each edge: 

= Y) / (djUj) ujdtj =-Y\( [djUjf dtj + Y][u j d j u j \\ 
Jo Jo 

3 3 3 

= -J2j o ( d 3 u jf dt 3 + ( d^^li ? Uj ^ ' 



independent of j * 



=0 



Clearly, $ £ Ai + (T), and any matrix belonging to Ai + (T) defines an eigenfunction to A = 0. Thus, the 
eigenspace Eo(T, A) is isomorphic to AA + (T) which shows the multiplicity formula for A = 0. 
Now suppose A > 0. In this case, a fundamental solution of (4.5) is given by 



(4.13) 
Thus, 



and, by (4.6), 



tt/ \ / /T^ sin(xv / A) T 
U(a;) = cos(a;V A)$ H I 7 . 

V A 

U'(a;) = -\/Asin(a;v / A)$ + cos(xVA)#, U'(l) = -VXsm(VX)$ + cos(v / A)* = 

sin \/A 



In the case 

we can conclude that 



U(l) = $* = $cos\/A- 
sin vA 7^ 0, 



v/A 



$ = 



1 



(cosVx * + J *A, 



\/A sin \/A ^ 

and, using $ e = and \f r = e* * A, we are led to 

(A* &/}*) e + cos\/A (yU^e*) e = 0, 

or the characteristic equation 

(4.14) 2(r)V> = - cos %/a v>- 

This shows the multiplicity formula for sin vA ^ 0. For the remaining case, suppose first 

cos V A = — 1. 

Then $* = — $, and (4.7 1 and (4.10 ) imply that m g (X) > N—n+1, using eigenfunctions with vanishing matrix 



ty. All eigensolutions vanishing in all vertices are of the form 



l(xy / A) 



* with a matrix * = vjr* = -0 e * * „4 



that has to be a multiple of A This shows m g (X) = N — n + 2. If instead 



cos 



Vx = i, 



then — <3>* = >]> = -0 e* * .4. can only be non-vanishing, and eigensolutions of the form sln ^ ^ can only 
exist if r is bipartite. On the other hand, $* = $, and (4.7) and (4.11) imply that the eigenspace E\(T, A) 
contains a subspace isomorphic to A^ + (r) using eigenfunctions with vanishing matrix "J. Thus, by ( |4.11 1, 
m g (X) = N — n + 2 in the bipartite case and m g (X) = N — n if T is not bipartite. □ 



Remark 4.2. 1) In the first part of Theorem ^. 1\ we cannot drop the assertion that X > 0, whether the graph 
is bipartite or not. This can be seen considering the following graphs 
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In the first case, — 1 is an eigenvalue of Z but is not an eigenvalue of A KC . In the second one, — 1 is 
not an eigenvalue of Z but the function of constant value (1,-1,0,-1,1) is an eigenfunction of A KC for the 
eigenvalue 0. 

2) Like in the case of CK -condition, the above result shows that in the relevant case of a connected graph 
the spectrum of A KC can be partitioned as 

<t(-A) = {OjU^Uff,, 

where 

a t := |(2£7r±arccosa) 2 | - a e a{Z) \ {-1,1} and I G zj , 
is the set of the immanent eigenvalues and the elements of 

<t s :={fcV|fceZ\{0}} 

correspond to the spectrum of the single edge problem. For the multiplicities of the eigenvalues we have 

(1) m(0) — N — n + 1 ifT is bipartite and m(0) = N — n ifT is not bipartite; 

(2) m(A) = dimker ^cos \/A — Z^j for A € c^; 

(3) m(fc 2 7r 2 ) = JV — n + 2 ifT is bipartite and m((2k + l) 2 7r 2 ) ~N — nifTis not bipartite; 

(4) m(4fc 2 7r 2 ) = N -n + 2. 

Remark 4.3. 1) Note that on trees there are no harmonic functions under the KC -condition. 

2) Bipartite graphs are exactly those graphs that can be endowed with an orientation such that each vertex is 
either a sink (only incoming incident edges) or a source (only outgoing incident edges). In this way harmonic 
functions with respect to the KC -condition are edgewise constant functions that satisfy the Kirchhoff rule at 
each node, i.e. C-flows in the sense of [321 §IX.4]. In fact, these edge distributions are exactly the elements 
of the null space of the incidence matrix T>(T). Now, it is well-known that the nullity of this matrix is 
N — n + c(r), cf. e.g. [5J §1.5]. This yields another proof for the claimed value of m g (0). 

3) In the non-bipartite case, the above characterization fails to hold since the Kirchhoff condition for the 
potentials on the edges does not take into account the orientation, while it is essential in the definition of the 
circuit space. To overcome this problem, observe that harmonic functions with respect to the KC-condition are 
edgewise constant functions that satisfy the Kirchhoff rule at each node ( regardless of the bipartiteness of the 
graph). In fact, their edge distributions are exactly the elements of the null space of the unsigned incidence 
matrix of V . Now, it is well-known that the nullity of this matrix is N — n + c + (T), cf. [32]. This yields 
another proof for the claimed value ofm g (0). 

In the same way as for the continuity ( |3.1[ ) under general consistent Kirchhoff flow conditions [4] , the above 
multiplicity formulae determine the asymptotic behavior of the eigenvalues of the elliptic problem (EPy^r). 

where R is a 27V x 2N hermitian matrix. 

Corollary 4.4. For the canonical Laplacian —A under the KC-condition the spectral asymptotics 

holds, where denotes the k-th eigenvalue, and where the eigenvalues are counted according to their (geo- 
metric = algebraic) multiplicities. 
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For a general, not necessarily hermitian matrix R, one can follow the Liouville transform approach in [6] 
counting the algebraic multiplicities in order to get the same asymptotic formula. We omit the details here. It 
has been shown in [2] that T is bipartite if and only if —1 is an eigenvalue of the matrix Z or, equivalently, the 
spectrum of Z is symmetric with respect to counting multiplicities. This can be used in various applications. 

Corollary 4.5. Suppose that the graph T is connected. Then the spectra of —A CK and —A KC coincide, 
counting multiplicities, if and only if T is unicyclic and bipartite. 



Proof. If T is unicyclic and bipartite, then corank(T) — 1 and the assertion follows by Theorems 3.2 and 4.1 
For the reverse implication, by hypothesis, the spectrum of Z is symmetric with respect to 0, which implies 
that r is bipartite. Moreover, since corank(r) = 1, i.e. N = n, and T contains exactly one circuit. □ 



In the disconnected case, both spectra coincide counting multiplicities, if and only if T is bipartite and each 
connected component is unicyclic, since iV = n if and only if corank(r) = c(F). 

Corollary 4.6. The graph T is bipartite if and only if the network immanent eigenvalues (sin^/X ^ 0) of 
-A CK and -A KC are the same, counted according to their multiplicities. 

Remark 4.7. IfT has a component with at least two edges, then Z induces an immanent eigenvalue between 
and the first non-trivial, non-immanent eigenvalue ir 2 . Since in the connected case the heat semigroup under 
the CK-condition always converges towards an equilibrium, while the semigroup under the KC '-condition 
only fails to do so if T is unicyclic with a circuit of odd length, this shows that the rate of convergence of 
the semigroup is usually determined by connectivity of the graph. Moreover, this eigenvalue coincides with 
/t := (arccosa) 2 and v := (arccos/3) 2 (for the CK-condition and the KC-condition, respectively), where 
a and (3 are the largest and the smallest eigenvalues of Z different from ±1, respectively. Observe that 
in the connected case, the eigenvalue 1 is necessarily simple by the Perron-Frobenius theorem. Hence, the 
eigenvalues of the dual pair (EPck), (EPkc) come in pairs and, in the particular case of bipartite graphs, 
even in quadruples (because then the spectrum of Z is symmetric with respect to the origin). 

As eigenvalues corresponding to ± cos are counted in both cases in the same order, the bipartite 

case permits the following comparison of the eigenvalues. For trees, we find 

Vfc G N : \k K < Af C - 



If r is bipartite and unicyclic, i.e. N = n with even circuit, then 



Vfc G N : Xk K = Af C 



If r is bipartite with general corank > 2 , then 

Vfc G N : X c k K > Af + V„ > Af c 
However, no general uniform comparison seems to be available in the non-bipartite case. 

Example 4.8. Let K n denote the complete graph on n > 3 vertices. Then the eigenvalues in the CK-case 
read 



1(5 
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Figure 1. On the abscissa, the eigenvalues of —A CK and —A KC are plotted in blue and 
red, respectively, in correspondence with the associated eigenvalues of Z on the ordinate axis, 
which are plotted in black. The first plot reflects the case of a non-bipartite graph, for which 
the spectrum of Z is not symmetric with respect to 0; while the latter corresponds to the 
bipartite case. 



(1) A = 


= 0, m(0) = 1, <p = e 




(2) A = 


= 47r 2 fc 2 , k ^ 0, m(A) 


= 2 + |n(n — 3), <p = e 


(3) A = 


:7 r 2 (2fc + l) 2 , m(A) = 


= \n{n — 3), <p = 


(4) cos 


v / A = -(n-l)" 1 , 


m(A) = n — 1, <£> G fcer e e* 



u>Me in i/ie KC-case we obtain the following ones. 
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(1)A = 


0, m(0) = |n(n-3), i> 


= 


(2) A = 


Ait 2 k 2 , k^O, m(A) = § 


ra(n — 3), ^ = 


(3) A = 


n 2 (2k + l) 2 , m(A) = 2- 


- |n(n — 3), ip = e 


(4) cos 


v/A = (?i- l)- 1 , m(A) 


= n — 1, ?/> € fcer e e* 



There is no uniform inequality between the X^ K and the Xj? c . For n — 3 e.g., Xq K < Xq , while X± > X^ . 

Remark 4.9. One of the most distinguished features of the characteristic equations for —A CK and —A K 
is the role played by Z. This becomes even more apparent once we recall that the spectrum of the transition 
matrix Z of a graph essentially agrees (up to a reflection in the point ~) with the spectrum of the so-called 
normalized Laplacian C defined by 

C := D-^LD-i = I — D~?AD~i = I - D^DZ D~i = (I - Z) D~* , 

where we have set 

L := Diag (A e) - A and D := Diag (A e) . 
Thus, the eigenvalues of Z are precisely those of I — £ counting multiplicities, in particular 
(4.16) A € cr(Z) «1-Ae er(£). 

This correspondence has been routinely exploited to study asymptotics of random walks on graphs, see e.g. [151 
§ 1.5], in particular because the matrix C has been studied very thoroughly over the last two decades. For our 
purposes, it is particularly interesting that applying known spectral comparison results for graph operations 
together with the strong monotony properties of the functions [0, it 2 } 3 X \— > ± cos \/A G [—1, 1], we can obtain 
some spectral monotony properties for —A CK and —A KC as well. 

Hence, some comparison results about the spectra of —A CK and —A KC become easily amenable, simply 
exploiting basic information about the spectra of C A curious examples is given in the following, where a 
special role is played by a strong kind of symmetry. We recall that a graph is said to be distance transitive if 
for any two pairs of vertices (x, y), (v, w) such that the distances between x, y and between v, w agree there is 
a graph automorphism mapping x to v and y to w, cf. 8, § 20]. E.g., complete graphs, cube graphs, and the 
Petersen graph have this property. 

Proposition 4.10. Assume V to be a distance transitive graph of diameter d. Then, both —A CK and —A K 
have exactly 2(d + 1) distinct eigenvalues in each interval {{2kir) 2 , (2(fc + l)ir) 2 }, k G N. 



Proof. Taking into account Theorems 3.2 and 4.1 the assertion is a direct consequence of [TH1 Thm. 7.10] (or 
of [SI Prop. 21.2], after observing that vertex transitive graphs are necessarily regular). □ 

Remark 4.11. If one allows Z and hence —A CK ,—A KC to be weighted, stronger assertions can be proved 
leading to spectral comparisons between Laplacians on distance transitive T and on some related path graphs, 
cf. |15[ Thms. 7.11 and 7.13]. We will not further elaborate on this point. 

Remark 4.12. Producing results analogous to that in Proposition 3.5 for —A KC is less obvious, since less 
information seems to be available about the smallest eigenvalue of Z, i.e., about the largest eigenvalue of C 
(apart from the bipartite case, of course). However, in the last few years many interlacing results have been 
obtained for the eigenvalues of £ under graph operations. For example, it is proved in [TU1 Thm. 1.2] that if 
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Ti is a graph on n vertices, 

T2 is a subgraph ofY\, 

t is the number of isolated vertices ofT 2 , 

T 2 is the complement ofT 2 in T\, i.e., the graph obtained by deleting from T\ the edges 0/T2, 



(1) J (2) 

i j ar^ " 
then for k = 1 



and [ij denote the eigenvalues of the normalized Laplacian ofTi, T 2 , respectively, 
,n 

(i) ( 2 ) J Mfc+t-i, */ r 2 is bipartite, 

I (J* k ! t , otherwise, 



where = . . . = /Aq := and /i^+i = • • ■ = Mn+t := 2- T^ 77 - analogous relation holds if a graph is 

added, instead of subtracted [101 Cor. 1.4] J. Yei more refined results related to more subtle structures (like 
coverings and spanning subgraphs) have been recently obtained in the doctoral thesis of Butler |12) and in still 



unpublished lecture notes by Hall [191 § 3]. In view of Theorems 3.2 and 4-1 it is clear how to translate all 



these interlacing results for the spectrum of C into interlacing results for spectral subsets of — A CK or —A KC . 

5. Inverse spectral aspects 

Just like under the CK -condition considered in Section[3] one cannot recover the graph from the eigenvalues 
of — A KC . In fact, one can use the same couple of regular isospectral, but non-isomorphic graphs as in [5] 
based on a well-known example by Hoffman and Ray-Chaudhuri, displayed by the left pair in Figj5] The 
right pair displays the same graphs where in each of both of them the white vertices have to be identified. 
They cannot be neither isometric nor isomorphic, since the left graph is planar, while the right one is not. 



Figure 2. Two non-isomorphic, but isospectral graphs under the K C-condition. 

Corollary 5.1. There are non-isomorphic pairs of regular graphs having the same point spectrum of the 
Laplacian —A under the KC -condition. 

In fact, the pair in Fig. [5] is not an exception. The above multiplicity formulae and the proof of Theorem 
|5.6| show also that the notions of isospectrality with respect the Laplacian under the CK and KC condition, 
which we denote by —A CK and —A KC throughout, are equivalent. 

Corollary 5.2. Two graphs Ti andY 2 are isospectral with respect to —A CK if and only if they are isospectral 
with respect to —A K . 

Remark 5.3. $From the eigenvalues of Z one can recover the number of vertices n and the bipartiteness, 
but, in general, the number of edges cannot be determined with the aid of o~(Z), see [3][TT]. A very simple 
pair of graphs that are isospectral with respect to Z, but differ in the number of edges is given by the circuit 



C4 and the claw ^1,3 (star with 3 edges) that have both the eigenvalues —2,0,0,2. By Theorems 3.2 and J^.l 
if two graphs are isospectral with respect to —A under the C K -condition or under the KC -condition, then 
they are also isospectral with respect to Z. But the latter condition is weaker than the first one, since one can 
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recover the number of edges from the —A-spectrum, see Theorem 5.6 Note that the pair C4, Ki^ is neither 
isospectral with respect to — A , nor with respect to —A KC , since the coranks are different. 

In the bipartite case, two vmicyclic graphs are isospectral with respect to —A CK if and only if they are 



isospectral with respect to — A KC , by Corollary 4.5 But, one cannot determine the graph by means of these 
eigenvalues in this class of graphs. 

Corollary 5.4. There are non-isomorphic pairs of bipartite unicyclic graphs that are isospectral with respect 
to —A under the C K -condition or under the KC -condition. 

Indeed, the following example by S. Butler and J. Grout [11] displays the desired properties. Let Ti be the 
circuit of length 8 and T 2 the circuit of length 4 with two 2-paths glued to one vertex as in Figure 3. Then 
the eigenvalues of both Z(Ti) and -Z(r 2 ) read 

1 — — — -1 

' 2 ' 2 ' ' ' 2 ' 2 ' 

The example also shows that one cannot recover the girth of the graph from the eigenvalues in question. 



Figure 3. The example from S. Butler and J. Grout [TT]: two non isomorphic bipartite and 
unicyclic, but isospectral graphs under the CK- and the if C-condition. 

Combining Corollaries |4.5| and |5.2| yields that there cannot be a pair of non-isomorphic and non-bipartite 
unicyclic graphs that are isospectral with respect to —A CK and —A KC . In other words 

Corollary 5.5. If two unicyclic graphs are isospectral with respect to —A CK or —A KC , and at least one of 
them is not bipartite, then they are isomorphic. 

Recall that the complexity k(T) of a graph T is defined by the number of spanning trees in T. By Kirchhoff 's 
matrix-tree theorem, it can be calculated from the eigenvalues of the combinatorial Laplacian L = D — A. 
But this is not possible from the spectrum of the matrix Z. In fact, a formula due to Runge and Sachs shows 
that knowing the determinant of Diag (A e) is necessary and sufficient in order to reconstruct the complexity 
of T from the spectrum of Z, cf. [TBI § 1-9, #10]. 

Theorem 5.6. From the eigenvalues (and their multiplicities) of either —A CK or —A KC , the following 
invariants can be determined: 

(1) the number of vertices n, 

(2) the number c(T) of connected components of the graph, 

(3) the number of edges N , 

(4) the numbers c + (r) and c~(T) of bipartite and non-bipartite components (and hence bipartiteness). 
However, it not possible to deduce from those eigenvalues the complexity, the degrees, and, in particular, 

whether the graph is regular. 

Proof. Under the Cif -condition, we first deduce c(T) = rn g (0) and 

c-(r) = i(m 9 (4^ 2 ) - m g (n 2 )) and c+(T) = c(T) - c"(r). 
Under the -fTC-condition, instead, we deduce 

c"(r) = ^K(7r 2 ) - m 3 (4^ 2 )) and c+(T) = m g (\ix 2 ) - m g (0), 
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and hence also c(T). Since Z is diagonalizable [2J, and since m g (l;Z) = c(T) and m g (—l;Z) = c + (T) (recall 
that r is bipartite if and only if — 1 g a(Z)), we deduce in both cases 

(5.1) n= m g (fx;Z) = c(T) + c + (r)+ £ m g (X). 

p,£<j(Z) Ae(0,ir 2 ) 



Using Theorems 3.2 and 4.1 we obtain N = n — 2c(T) + m g (Air 2 ) in the first case, and N — m g (Q) + n — c + (T) 
in the second one. 

The pair of isospectral, but non-isomorphic graphs in Figure 3 shows that it is impossible to determine 
the vertex degrees, in particular whether they are all equal, from the data in question. Finally, the same pair 
displays complexities n = 8 and k = 4, respectively. □ 

Remark 5.7. It is well-known that n, N , and c(T) can be deduced from the adjacency matrix A, or from a 
combination of both the combinatorial and the normalized Laplacian. Moreover, it is known that c + (T), the 
number of bipartite components, agrees with the multiplicity of as an eigenvalue of the signless Laplacian 

Q := Diag (A e) + A, 

a matrix which has gained much attention in recent years, cf. [17] . To the best of our knowledge, however, it 
is unknown whether c + (T) can be found only knowing the spectrum of any of the other relevant matrices; and 
also whether c(T) can be reconstructed from the knowledge of the spectrum of Q. Hence, it seems that A CK 
and A KC are, taken individually, more comprehensive sources of spectral information than the usual matrices 
considered in spectral graph theory. 

In the regular case (say, with degree 7), the notions of isospectrality with respect to 

• the adjacency matrix A, 

• the combinatorial Laplacian L = D — A = jl — A, 

• the transition matrix Z := D~ X A = r y~ 1 A, 

• the normalized (discrete) Laplacian C :— D~zLD~2 = I — j~ 1 A, and 

• the signless Laplacian Q := D + A = 7/ + A, 

where D := Diag (A e), are mutually equivalent. Thus, we are led to the following 



Corollary 5.8. IfY is regular, then one can recover all the invariants mentioned in Theorem \5.6[ as well as 
the complexity from the eigenvalues of each of the operators or matrices —A KC , —A CK , A, L, Z, C, and Q. 



Remark 5.9. Combining the multiplicities for harmonic functions, Theorems ^. 1\ and \3.£\ yield 

dim Ker A KC - dim Ker A CK = N — n — c~(T). 

In the cohomological language, —A KC can be written as dd* and —A CK as d*d, where d denotes the first 
derivative with continuity vertex conditions. It follows that 

ind d := dim Ker A KC - dim Ker A CK 

agrees with the Euler characteristic (N — n) ofT ifT is bipartite (c~ (T) = 0). This yields an easy vroof of [2U1 
Thm. 20] in the bipartite case, but only then. The formula ind d = N—n does not hold in the non-bipartite case. 
This should be compared with the theory developed in |30j . The result of Fulling, Kuchment and Wilson aims 
at comparing the null spaces of an operator acting on L 2 (0,1; £ 2 (E)) -functions and another operator acting 
on the space of their derivatives (i.e., on a space of 0-forms and 1-forms, respectively), whereas we regard both 
operators as acting on the same space. If the underlying graph is bipartite, then it is always possible to find 
an orientation that allows to formulate KC-condition for a function u as the CK-condition for its derivative, 
as done in [20]; in the non-bipartite case, this is in general impossible, as even the example of a circuit of 
length 3 shows. This explains why our results are not comparable with those of [20] in the non-bipartite case. 
In particular, their formula (36) and our Theorem 5. 1 are only seemingly mutually contradictory - in fact, 
they refer to different objects. 
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